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Abstr act. Wc investigate quantum control of the dissipation of entanglement 
under environmental decoherence. We show by means of a simple two-qubit 
model that standard control methods - coherent or open-loop control - will not 
in general prevent entanglement loss. However, we propose a control method 
utilising a Wiseman-Milburn feedback/measurement control scheme which will 
effectively negate environmental entanglement dissipation. 
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1. Introduction 

Entanglement has recently emerged as a significant resource in quantum infor- 
mation and other applications. However, an important problem is to ensure the 
robustness of this resource; that is, the ability to maintain it against decay. Un- 
like other quantum properties, entanglement is not invariant under general unitary 
transformations, and this implies the possibility of decay under unitary forces. 
These arc however, reversible and thus conversely enable entanglement production, 
given appropriate control methods. 
^-H Additionally, environmental dissipation is an ever-present source of the loss of 

entanglement; and it is therefore important to devise quantum control procedures 
to protect against this loss, where possible. 

In this note we analyze the possibility of preservation of entanglement against 
''"^ decay for a simple two-qubit Bell state by means of quantum control. We con- 

elude that this is not possible for a simple hamiltonian quantum control process 
>> (open-loop control) but show that an appropriate feedback/measurement control 

k>( procedure can prove effective protection against environmental loss in certain cases. 
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2. Unitary Dissipation 

Quantum dissipation is usually associated with non-unitary, non-reversible pro- 
cesses. However entanglement is subject to unitary dissipation, since unitary evo- 
lution associated with a (hermitian) hamiltonian does not necessarily preserve en- 
tanglement. Conversely, entanglement may be produced by the evolution induced 
by a quantum control hamiltonian. We choose a simple example to illustrate the 
bipartite case. 
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2.1. Entanglement production. Consider the unitary evolution U{t) induced by 
the hamiltonian H given by 



(2.1) 
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which corresponds to a system with Heiscnberg coupling and local control terms 

(2.2) H = aZ(E)I + bI(E)Z + c{X(x)X + Y(x)Y + Z(E)Z) 

with xi — a + b + c, X2 = a — b — c, = —a + b — c, x^ = —a — b + c, y = 2c. For 
simplicity we assume a ~ b so that X2 ~ x^. 

We shall use the Concurrence C jlj as a measure of entanglement for a bipartite 
two-qubit system. The concurrence C of a two-qubit state p is given by 

(2.3) C max{Ai - A2 - A3 - A4,0}, 

where the quantities Ai are the square roots of the eigenvalues of the 4x4 matrix 

(2.4) p{Y(g)Y)p*{Y(g)Y) 

"0 - 



in descending order, where Y 







Acting with the unitary evolution matrix U{t) — exp(itH) induced by Eq.(2.1) 
on the base vector vq = [0, 1, 0, 0]^ gives 

v(t) — exp {itH)vQ 

= e**"^ [O,cos{ty),ism{ty),0f, 



from which we can easily obtain the concurrence of v{t) to be | sui{2ty)\. 



2.2. Unitary dissipation of Entanglement. The off-diagonal coupling term y 
changes the entanglement. This can be used to create entanglement but it also 
destroys entanglement. 

Referring to Figure 1, we see immediately that if we start in the state vq then at 
t = 7r/4 (in units ofl/y) we have the maximally entangled (Bell) state -^[0,1, i,0]^ 
(up to a global phase factor) but the unitary action U{t) destroys the entanglement, 
completely a,tt = tt/2. Similarly, if the system is prepared in a maximally entangled 
Bell state such as 

(2.5) |*Bcii) = ;^[0,1,1,0]^ 

then the entanglement is not preserved unless y — 0, which corresponds to no 
coupling. 



Note that for such calculations it is important to choose a fixed basis - here we choose the 
standard basis. 
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Figure 1. Concurrence vs Time t (units of l/y) 



3. Environmental Dissipation 

The standard form of a general Markovian dissipative process in Quantum Me- 
chanics is governed by the Liouville equation obtained by adding a dissipation 
(super-)operator LD[p(t)] to the usual hamiltonian term: 



(3.1) 



p{t) = -^[H,p{t)]+LD[p{t)] 



where the density matrix p for an iV-level system is an x (semi-)positive 
matrix. 



3.1. Lindblad Equation. In order that the quantum Liouville equation Eq.(3.1) 



should define a physical dissipation process, certain constraints must be imposed, 
most notably complete positivity. The appropriate constraints emerge from physi- 
cal stochastic dissipation equations such as those given by Lindblad and others in 
differential form ^ , as well as in global form [5] . The result is that completely pos- 
itive evolution of the system is guaranteed provided the dissipation super-operator 
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Ld takes following form 

(3.2) LuW)] = - Y^{[V,p{t), + [F„ p{t)V}]} 

s=l 

where the matrices Vg are arbitrary. Denoting the standard basis for NxN matrices 

by 

(3.3) {Eij)mn = SimSjn {i, j , ui, n = 1 . . . N) 

and relabelling, using the notation [m; n] = (m — 1)N + n, we have 

(3.4) V[i.j] ^ a[,.j]E[,.j]. 



The virtue of Eq. (3.1) is that every Markovian dissipation process has to satisfy 
it and so results derived from its use have great generality. A drawback, how- 
ever, is that it is often not obvious how to relate phenomcnological observations of 



dissipative effects to dissipation operators of the form (3.2) 



3.2. Phenomenological description. Having chosen a certain computational or 
preferred basis, one can phenomenologically distinguish two types of dissipation: 
phase decoherence and population relaxation. The former occurs when the interac- 
tion with the environment destroys the phase correlations of certain superposition 
states. In the simplest case this leads to a decay of the diagonal elements Pkn{t) of 
the density operator at a constant (dephasing) rate Tkn'- 

(3.5) Pkn{t) = -i{[H, p{t)])kn - ^knPkn it). 

Population relaxation occurs, for instance, when a quantum particle in a certain 
state spontaneously emits a photon and transitions to a less energetic quantum 
state. In the simplest case, when there are only jumps between the basis states, \k) 
and \n) say, occurring at fixed rates jkn, the resulting population changes can be 
modelled as 

(3.6) pnn{t) = -l{[H, p{t)])nn + 

nk Pkk{t) 

k^n 

where jknPnn is the population loss for level \n) due to transitions |7i) —> and 
InkPkk is the population gain caused by transitions \k) — >■ \n). The population 
relaxation rate ^kn is determined by the lifetime of the state \n) and, for multiple 
decay pathways, the relative probability for the transition \n) — |fc). 

Phase decoherence and population relaxation lead to a dissipation super-operator 
(represented by an iV^ x N'^ matrix) whose non-zero elements are 

(3.7a) {LD)[k-n]\k-n\ = -rfen fc ^ 71 

(3.7b) {LD)[n-n],[k-k\ = +7«fc k^n 

(3.7c) {LD)[n-n],[n-n] = - ^ Ikn 

n^k 

where Tkn and jkn are taken to be positive numbers, with Tkn symmetric in its 
indices, and again we employ the convenient notation [m;n] = (m — l)iV + n 
introduce above. 

The N"^ X matrix super-operator Ljj may be thought of as acting on the 
A''^-vector r obtained from p by 

(3.8) ^[m;n] = Pmn- 
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The resulting vector equation is 

(3.9) i = Lv = {Lh + Ld)v 

where Lh is the anti-hermitian matrix corresponding to the hamiltonian H . We 
obtain Lh explicitly by using the standard algebraic trick applied in evaluating 
Liouville equations (see, for example [2]). The correspondence between p and r as 
given in Eq. (3.8 1 tells us, after some manipulation of indices, that 

(3.10) p^v^ ApB^ A®Bv 
using the direct (Kronecker) product of matrices. 

3.3. Constraints on Dephasing Rates. The phcnomenological description above 
does not impose any constraints on the population relaxation and decoherence pa- 
rameters present in the dissipation matrix. In practice, however, the values of the 
dissipation parameters F^n and ^kn must satisfy various constraints [3] to ensure 
that they describe physical processes, which can be derived by use of the Lindblad 
equation. 

For simplicity we restrict ourselves here to the case of pure dephasing. Ex- 
perimentally, this is often the dominant decoherence process as the population 
relaxation (or Ti) times for most systems are much longer than the dephasing (or 
T2) times so that we may effectively neglect the relaxation rates 7. In the pure 
decoherence (dephasing) case, comparison of Eq.(3.2) and Eq.(3.4) with Eq.(3.7) 
tells us that the 7 terms vanish if we choose aj^.j] = for i 7^ j. The decoherence 
parameters Tij are then given by 

(3.11) r,, = i(|a[,;,]p + |a[,;,]n {i,j^l...N i^j). 

This leads to a mathematically very simple situation, as the dissipation matrix 
Ld{) is then diagonal. For the = 4 system, this gives six pure dephasing param- 
eters iVij = Fji, Tii — 0), determined by four constants, so there are two relations 
between the F's. Explicitly, we have for the two-qubit, 4-level case, 

Ldo diag{0, -F12, -F13, -ri4, -Fai, 0, -Fas, -F24, 

(3.12a) - F31, -F32, 0, -F34, -F41, -F42, -F43, 0} 

with Tij = Tji, and the constraints that must be imposed to ensure that we have 
physical process are then 

(3.13) F12 -f F34 = Fi4 -I- F23 = F13 + F24. 

3.4. Markovian dissipation of entanglement. Wc now give an example of a 
standard dephasing process acting on a maximally entangled state. 

Example 1. Consider the Bell state 

(3.14) i;b = ^[0,1,1,0]^. 

The corresponding Liouville vector is r = \ [0, 0, 0, 0, 0, 1 , 1, 0, 0, 1, 1, 0, 0, 0, 0, 0] 
and the action of the dephasing operator L/jo of Eq.(3.12| is given by 

(3.15) r = ir = (L,3o)r, 



T 



as in Eq. (3.9). This equation may be immediately integrated to give 
1, 



(3.16) v{t) = -[0,0,0,0,0, l,e~^^^*, 0,0,6-^^3*, 1,0, 0,0, 0,0] 
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concurrence 




Figure 2. Concurrence vs Time t (units of I/T23) 



corresponding to the density matrix 



(3.17) 
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Note that this does not represent a pure state except at t = 0. The concurrence as 



defined in Section 2.1 evaluates to exp(— r23t). (See Figure 2). 



The results of Example [T] are essentially unchanged in the presence of an ad- 
ditional free (i.e. diagonal) hamiltonian as this commutes with the dissipation 
super-operator LdoE]- 
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4. Quantum control of entanglement dissipation 



In this section wc analyze two quantum control schemes for mitigating entangle- 
ment dissipation. Wc first consider a simple scheme based on open-loop coherent 
control and then a measurement-based feedback scheme. 

4.1. Coherent quantum control. For pure dephasing and a hamiltonian of the 
form (2.11 the (2, 3)-subspace is invariant, i.e. if a state starts in this subspace it 



remains there. Thus if the system is initially prepared in the Bell state Eq. (3.14) 
then we may consider the reduced dynamics on the (2, 3)-subspace. The density 
operator p restricted to the (2, 3)-subspace is essentially a qubit, i.e., a 2 x 2 positive 
matrix p. Introducing the Pauli matrices 



(4.1) 



X = 



1 

1 



Y = 







Z = 





-1 



then the Lindblad equation Eq(3.2|, with a single V term — Z, takes the simple 
form 



(4.2) 



f>^-i[yX, p]+^{Z)p 



where D{Z)p = ZpZ - p. 

It is convenient to use the Bloch representation, for which the 3-vector s associ- 
ated with the density matrix p is given by 



(4.3) 



(trace(Xp), trace(yp), trace(Z/})) 



For hamiltonian dynamics, ||s|| is constant and the motion is governed by the or- 
thogonal group 0(3) acting on the Bloch Sphere. In the case of dissipation, the 
motion takes place in its interior, the Bloch Ball, and is locally governed by the 
afhne Lie algebra (?/(3) ® R^; this is globally a semi-group due to boundary condi- 
tions \J\. Therefore in general the Bloch equation is 

(4.4) s{t) = As{t) 



c. 



In the simple case considered here, we have 



(4.5) 



A = 



-2 
-2 
2y 








For 1/ 7^ the Bloch matrix A is invertible and any initial state eventually goes to the 



unique steady state Sg 



-A c = [0,0,0] , which corresponds to the completely 



mixed state on the (2, 3)-subspace, which has zero concurrence. Thus regardless of 
the value of the control y the system will eventually go to the completely mixed state 
on the (2, 3)-subspace and all entanglement will be lost. If y = (no control) then 
all the Bloch vectors along the z-axis in the (2, 3)-subspace are steady states, but if 
we start with a Bell state we still go to the completely mixed state on the subspace 
for i — > oo. So in this case there is no way coherent control can prevent the decay of 
the entanglement. This is understandable as the hamiltonian term determines the 
anti-symmetric part of A while the decoherence term defines the symmetric part, 
so the control cannot change the contraction of the Bloch vector introduced by the 
dephasing term. If we want to stabilize an entangled state we need to employ a 
more sophisticated method. One approach is to utilize a measurement and feedback 
scheme. 
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Figure 3. Basic direct-feedback control setup. 



4.2. Measurement /feedback control scheme. One approach that has shown 
recent promise in particular for stabilizing quantum states (8] is reservoir engineer- 
ing, in particular using direct feedback. 

Suppose we have a system whose evolution is governed by a Lindblad master 
equation (3.1). If we add a fixed continuous weak measurement of an observable 



M and apply a fixed feedback hamiltonian F conditioned directly on the measure- 
ment current, then according to the general theory developed by Wiseman and 
Milburn [3J [TO] the master equation is modified as follows: 



(4.6) 



where H is the original system hamiltonian plus the open-loop control hamiltonian 
i?o + tic plus a feedback correction term {AI^F + FM)/2 (see the diagrammatic 
scheme of Figure 3) . 

Suppose we have a system with hamiltonian _ffo as in (2.1) subject to environ- 
mental dephasing described by a diagonal Lindblad operator V . If we now add a 
weak continuous measurement of the observable ^/mZ ® I and apply a feedback 
hamiltonian of the form ^/fX (E) X then the evolution of the system according to 



(4.6) is governed by 



(4.7) 



-i[H,p]++D[V]{p)+D[M ~iF]p 



As observed before, the (2, 3) subspace is invariant under H and V. It is also 
invariant under M — iF, i.e. any state starting in this subspace will remain in it. 
Thus, if we are only interested in initial states in this subspace we can again restrict 
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our attention to the dynamics on this subspace. The subspacc operators are 

(4.8) H^„^-''^ ^ nZ + yX, M^^'^^ = V^Z, F^^-S) = //X, V'^^^''>=VtZ 

and noting that Z^X + XZ = shows that the master equation for the density 
operator p restricted to the (2, 3) subspace takes the form 

(4.9) ^ = -i[yX + iiZ,p] +T>[V^Z - i^lA'P + rS)[Z]p. 

Here p, is the effective energy level splitting and F the effective environmental doco- 
herence rate, and y, yjm and -^7 are the effective strengths of the open-loop control, 
measurement and feedback, respectively. The corresponding Bloch operator is 

/m + r p 0> 

(4.10) A = -2\ -fj. f + m + T y 

V -y fj 

and c = [0,A^/mf,0\'^. 

For calculational simplicity we choose the open- loop control term to be 0, i.e., 
y = 0. In this case it is easy to compute the eigenvalues of A: 

(4.11) -2/,-/-2r-2m± V/'- V- 

Thus A is non-dcgcncratc if / 7^ and the system therefore has a unique steady 
state Sss = —A~^c on the subspace, which in density operator form is 



(4.12) p,, = 



1 ^/7^{iJ.+i{r+m)) 

2 ^2+(r+TO)(r+m4-/) 
v7^(M-»(r-i-m)) 1 

At2 + (r+m)(r+m+/) 2 



The purity of the steady state is given by 

(413) P -1+2 W + + 

^ ' ~ 2 + ^ (M^ + (r + m)(r + m + f)^ 

and the concurrence is 

UU) c o V^(M^ + (r + m)^)^/^ ^ 

^^•^^^ ^^^-%2 + (r + ^)(r + m + /) 

Note that C^s = yf^Pss — Ij i-e. there is a one-to-one correspondence between the 
concurrence and purity in this case. If we further choose = (adjust energy level 
splitting to be 0) then the expressions simplify: 

/to 



(4.15) - ' ' 



\~*r+/+TO 

and the purity and concurrence are given by 

1 2/m 
(4.16) P„ ■ 

(4.17) 




2 (r + / + m)2 



r + m + / 

For a given fixed dccohercncc rate F > 0, we can choose to and / in units of F 
and set F = 1. In this case the steady-state concurrence becomes 



10 



SOPHIE G. SCHIRMER AND ALLAN L SOLOMON 



log^Q(1-2(mf)i"^/(1+m+f)) 




10 20 30 40 50 60 70 80 90 100 



Figure 4. Steady-state concurrence Css, or precisely log]^Q(l — 
Css), as a function of the measurement and feedback strength m 
and /, respectively, assuming F = 1. 

Figure |4] shows the logarithm of 1 — Css as a function of the measurement and 
feedback strengths. We see that the optimum choice to maximize the concurrence 
is m — f with m and / as large as practically feasible. The plot given in Figure 
4 shows that if the measurement and feedback strengths are about 100 times the 
environmental decoherence rate F then the steady-state concurrence is greater than 

1-10-2. 

Although this scheme may be difficult to realize in practice because we require 
a non-local feedback hamiltonian of the form X <Si X, with some experimental 
ingenuity such a procedure may well be implemented. 
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